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I. INTRODUCTION. LEADING LOGARITHMIC CONTRIBUTIONS OF ORDER 

a'^{Za){m/M)EF 

It is well known that the radiative-recoil corrections of order a'^{Za){m/M)E^ to hy- 
perfine splitting in muonium are enhanced by the large logarithm of the electron-muon mass 
ratio cubed 0. The leading logarithm cube contribution is generated by the graphs in Fig. |l| 
with insertions of the electron one-loop polarization operators in the two-photon exchange 
graphs. It may be obtained almost without any calculations by substituting the effective 
charge a{M) in the leading recoil correction of order {Za){m/M)Ep, and expanding the 
resulting expression in the power series over a M. 
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FIG. 1. Graphs with two one- loop polarization insertions 



FIG. 2. Graphs with two-loop polarization insertions 

Calculation of the logarithm squared term of order a'^{Za){m/M)Ep is more challenging 
. All graphs in Figs. |l|, ||, H, ^, and |^ generate corrections of this order. The contribution 
induced by the irreducible two-loop vacuum polarization in Fig. § is again given by the 
effective charge expression. Subleading logarithm squared terms generated by the one-loop 
polarization insertions in Fig. |l| may easily be calculated with the help of the two leading 
asymptotic terms in the polarization operator expansion and the skeleton integral. The 



^We define the Fermi energy as 

E,='izv!^r^ycHB^, (1) 

S M \m J 

where ra and M are the electron and muon masses, a is the fine structure constant, c is the velocity 
of light, h is the Planck constant, i?oo is the Rydberg constant, and Z is the nucleus charge in 
terms of the electron charge [Z = 1 for muonium). The Fermi energy Ep does not include the 
muon anomalous magnetic moment a^ which does not factorize in the case of recoil corrections, 
and should be considered on the same grounds as other corrections to hyperfine splitting. 

^And by the diagrams with the crossed exchanged photon lines. Such diagrams with the crossed 
exhanged photon lines are also often omitted in other figures below. 
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FIG. 3. Graphs with radiative photon insertions 
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FIG. 4. Graphs with polarization insertions in the radiative photon 

logarithm squared contribution generated by the diagrams in Fig. ^ is obtained from the 
leading single-logarithmic contribution of the diagrams without polarization insertions by 
the effective charge substitution. An interesting effect takes place in calculation of the 
logarithm squared term generated by the polarization insertions in the radiative photon in 
Fig. |. One might expect that the high energy asymptote of the electron factor with the 
polarization insertion is given by the product of the leading constant term of the electron 
factor — 5a/(47r) and the leading polarization operator term. However, this expectation 
turns out to be wrong. One may check explicitly that instead of the naive factor above 
one has to multiply the polarization operator by the factor — 3a/(47r). The reason for this 
effect may easily be understood. The factor — 3a/(47r) is the asymptote of the electron 
factor in massless QED and it gives a contribution to the logarithmic asymptotics only after 
the polarization operator insertion. This means that in massive QED the part — 2a/(47r) 
of the constant electron factor originates from the integration region where the integration 
momentum is of order of the electron mass. Naturally this integration region does not give 
any contribution to the logarithmic asymptotics of the radiatively corrected electron factor. 
The least trivial logarithm squared contribution is generated by the three-loop diagrams 
in Fig. ^ with the insertions of the light by light scattering block. Their contribution was 
calculated explicitly in p[. Later it was realized that these contributions are intimately 
connected with the well known anomalous renormalization of the axial current in QED [Q. 
Due to the projection on the HFS spin structure in the logarithmic integration region the 
heavy particle propagator effectively shrinks to an axial current vertex, and in this situation 
calculation of the respective contribution to HFS reduces to substitution of the well known 
two-loop axial renormalization factor in Fig. ^ [Q in the recoil skeleton diagram. Of course, 
this calculation reproduces the same contribution as obtained by direct calculation of the 
diagrams with light by light scattering expressions. From the theoretical point of view it is 
interesting that one can measure anomalous two-loop renormalization of the axial current 
in the atomic physics experiment. 

The sum of all logarithm cubed and logarithm squared contributions of order 
a^{Za){m/M)EF is given by the expression [IIIH 
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FIG. 5. Graphs with light by light scattering insertions 



FIG. 6. Renormalization of the fifth current 

,^ / 4,nM A^ ^M\ a^(Za)m ~ 
\3 m 6 m / TT-^ M 

It was also shown in that there are no other contributions with the large logarithm of 
the mass ratio squared accompanied by the factor a^, even if the factor Z enters in another 
manner than in the equation above. 

Single-logarithmic and nonlogarithmic terms of order a'^{Za){m,/M)EF are generated 
by all diagrams in Figs. |I]-^ by the graphs with the muon polarization loops, by the graphs 
with polarization and radiative photon insertions in the muon line, and also by the graphs 
with two radiative photons in the electron and/or muon lines. Only a partial result for the 
single-logarithmic and nonlogarithmic corrections generated by the pole part of the graphs 
with both electron and muon polarization loops is known now |^. Numerically respective 
contribution is about 9 Hz, and may be considered only as an indication of the scale of 
the respective corrections. Corrections of this scale are phenomenologically relevant for 
modern experiment and theory 0. In this paper we calculate all radiative- recoil corrections 
generated by the diagrams including only the polarization loops, either electronic or muonic, 
leaving calculation of the other contributions for the future. 



II. TWO-PHOTON EXCHANGE DIAGRAMS. CANCELLATION OF THE 

ELECTRON AND MUON LOOPS 

Calculation of single-logarithmic and nonlogarithmic radiative-recoil corrections of rela- 
tive order a'^{Za){m/M) (and also of orders {Z'^aY{Za){m/M) and a{Z'^a){Za){m/M)) 
resembles in many respects calculation of the corrections of relative orders a{Za){m/M) 
and Z'^a{Za){m/M). It was first discovered in [0,|| that the contributions of the diagrams 
with insertions of the electron and muon polarization loops partially cancel, and, hence, it 
is convenient to treat such diagrams simultaneously^. Similar cancellation holds also for the 
corrections of order a'^{Za){m/M)EF, so we will first remind the reader how it arises when 



^We always consider the external muon as a particle with charge Ze, this makes origin of different 
contributions more transparent. However, somewhat inconsequently we omit the factor Z in the 



one calculates the polarization contribution of order a{Za){m/M)Ep. The nonrecoil contri- 
bution in the heavy particle pole of the two-photon exchange diagrams exactly cancels in the 
sum of the electron and muon polarizations (see for more details [§,§). Then the skeleton 
recoil contribution to the hyperfine splitting generated by the diagrams with two-photon 
exchanges in Fig. |^ is the result of the subtraction of the heavy pole contribution 
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FIG. 7. Diagrams with two-photon exchanges 
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where fj, = m/(2M), and 
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The electron polarization contribution is obtained from the skeleton integral by multi- 
plying the expression in eq.(0) by the multiplicity factor 2, and inserting the polarization 
operator {a/Ti)k'^Ii{k) in the integrand 
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The muon polarization contribution is given by a similar expression, the only difference is 
that 



I,{k) ^ h,{k:) ^ f 
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(6) 



Then the total recoil contribution induced by the diagrams with both the one-loop electron 
and muon polarizations in Fig. ^ has the form 



^^ = ^^^M^'Jo T 



fif^k) - f 



[eh{k) + eh^{k)]. 



(7) 



Next we rescale the integration variable k — > kM/m in the muon term and obtain 



case of the closed muon loops. The reason for this apparent inconsistency is just the cancellation 
which we discuss now. 
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FIG. 8. Diagrams with one-loop polarization insertions 
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We see that the electron and muon polarization contributions have partially cancelled. More- 
over, it is not difficult to check explicitly that the term with f{k/{4fi)) generates only cor- 
rections of higher order in /x, so with linear accuracy in the small mass ratio m/M all recoil 
contributions generated by the diagrams with the one-loop electron and muon polarization 
insertions in Fig. |^ are given by the integral 



^ ^ a(Za) m ~ 
AE = 8 ^ , ^ —Ep 



TX" 



M 



Hk 
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(9) 



This integral was calculated in |§] and we will not discuss its calculation here. Our only goal 
in this Section was to demonstrate the mechanism of the partial cancellation of the electron 
loop and muon loop contributions. 

III. DIAGRAMS WITH EITHER TWO ELECTRON OR TWO MUON LOOPS 

The nonrecoil contribution generated by the diagrams with two electron or muon loops in 
Fig. ID and Fig. ^ was obtained a long time ago ^. Although it was not emphasized in that 
work explicitly, it is easy to check that the result in [^ includes heavy pole contributions 
which are due to the diagrams with both the electron and muon polarizations. Repeating the 
same steps as in the previous Section, it is easy to see that the recoil contribution generated 
by the diagrams in Fig. |l] and Fig. |^ is determined by the integral 
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FIG. 9. Graphs with two muon one-loop polarization insertions 
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(10) 



where the numerical factor before the integral is due to the multiplicity of the diagrams, 
and the whole integral is similar to the integral in eq.(P). The only significant difference is 



6 



that now we have the two- loop factor k^P{k) in the integrand instead of the one-loop factor 
eh{k). 

We calculate the integral in eq. ( [T0| ) separating the contributions of small and large mo- 
menta with the help of the auxiliary parameter a such that 1 -C cr -C l//i 






Ep. 



(11) 



Then for the small integration momenta region in the leading order in /icr {^k < ^cr <^ 1) 
we have 



'^ dk _ 



^^' =^[ f /(/^^)^'A^(^) - -3/^'^ -^kHl{k). 



(12) 



We substitute in this integral the closed expression for the polarization Ii{k), and again 
preserving only the leading contributions in yucr ^ 1 obtain 
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(13) 



The high momenta contribution is calculated by expanding the polarization operator in 
1/k^ < l/a^ < 1 
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For calculation of this integral we use the standard integrals introduced in [|I2] as well as 
some new standard integrals (see Appendix), and obtain 
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Now we are ready to write down the total recoil contribution generated by the diagrams in 
Fig. I and Fig. | 
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The logarithm cube and logarithm squared terms in this expression are already known [|T],3 , 
and the single-logarithmic and nonlogarithmic terms are obtained here. 



FIG. 10. Graphs with both the electron and muon loops 



IV. DIAGRAMS WITH BOTH THE ELECTRON AND MUON LOOPS 



Consider now the diagrams with one electron and one muon loop in Fig. ^. We can look 
at these diagrams as a result of the electron polarization operator insertions in the muon 
loop diagrams in Fig. ^ The complete analytic expression for the last two diagrams in Fig. 
Bl has the form 
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where 



f{k) = f{k) + 



(17) 



(18) 



Unlike eq.(0) we have restored in eq.(p!7D the heavy particle pole contribution, which in 
the case of the muon polarization loop also carries the recoil factor. To simplify further 
calculations we rescale the integration momentum k -^ kM/m 
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and note that with the linear accuracy in m/M we may omit the second term in the square 
brackets in the integrand. Then the muon loop diagrams in Fig. ^ are described by the 
expression 



A^ = 8 



a{Za) m ~ ['^ dk~fk\ ^ 



(20) 



The integral in eq.fPOD turns into the contribution of the diagrams in Fig. [1^ after multipli- 
cation by the factor 3 and insertion in the integrand of the additional factor 
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This extra factor enters in the asymptotic regime since the characteristic scale of the inte- 
gration momenta in eq. f^Of ) is about one, and the parameter /x goes to zero. 

Then the contribution to HFS of the diagrams in Fig. [TO is given by the integral 
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After calculation we obtain (see Appendix 

AE = 
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V. DIAGRAMS WITH SECOND ORDER POLARIZATION INSERTIONS 

The recoil contribution to HFS generated by the diagrams in Fig. ^ and Fig. |ll] with 
two- loop electron and muon polarization insertions is given by the integral (compare eq.(P)) 



FIG. 11. Graphs with muon two- loop polarization insertions 
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where {a^ / Ti'^)k'^ l2{k) is the two-loop polarization operator |TT,T2 
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To simplify further calculations we represent the two-loop polarization operator in the 
form 



I^{k) = lh{k)^ f'dv 
4 Jo 
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where 
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The integral in eq.(^) decreases as 1/fc^ at large k. Note that since -R(f) -^ 3(1 — v) at 
V —^ 1 this leading term in the asymptotic expansion is not enhanced by the large logarithm 
In k. Absence of this logarithm significantly simplifies further calculations. 

In terms of the function R[v) the integral for the recoil contribution in eq. (p^ has the 
form 
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The first contribution on the right hand side is proportional to the well known one-loop 
contribution in eq.(|) |3|,|3,i 
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Calculation of the second term AE'^ is a bit more involved. We again introduce the 
auxiliary parameter a (1 ^ a ^ Va*) ^^'^ consider separately the small and large momenta 
contributions. For the small integration momenta region in the leading order in fia we have 
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where we used certain integrals for the function R{v) collected in the Appendix. 
The high-momentum contribution is given by the integral 
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First we use that A;^ ^ a^ ^ 1 and that R{v) — * 3(1 — f) asf — *> 1 and omit 4 in the 
denominator in the integrand, and then perform the calculations using the integrals from 
the Appendix 
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The total recoil contribution to HFS generated by the diagrams with two-loop polariza- 
tion insertions in Fig. |ll] is given by the sum of the contributions in eq.(^), eq.(|30D, and 
eq.m 



AE = AE" + AE 
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The logarithm squared term in this expression was obtained in 0, and the single-logarithmic 
and nonlogarithmic terms are obtained here. 

VI. DISCUSSION OF RESULTS 



Collecting all contributions in eq. (pIBf) , eq.(^), and eq.(|33D we obtain 
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The factor in the square brackets is about (—103), and numerically the respective con- 
tribution to the muonium HFS is 

AEnew = -0.027 7 kHz. (36) 

The magnitude of this correction is just in the range we should expect based on the 
partial result in p. The contribution in eq.(|36|) is of the same scale as the logarithmic 
in Za corrections of order {Za)^{m/M)Ep and a{Za)'^{m/M)Ep, calculated recently in 



14,15 



Collecting the recent results from 1 10, 14, 15] and eq. (p5D , and using the experimental 



value of the muonium hyperfine splitting [Q we may derive a value of the electron-muon 
mass ratio 

— = 206.768 279 8 (23) (16) (32), (37) 

m 

where the first error comes from the experimental error of the hyperfine splitting measure- 
ment, the second comes from the error in the value of the fine structure constant a, and the 
third is an estimate of the yet unknown theoretical contributions. 
Combining all errors we obtain the mass ratio 

— = 206.768 279 8 (43), 6 = 2- 10-^ (38) 

m 



which is almost six times more accurate than the best direct experimental value in |jT6|. 



Estimating the errors in eq.(p7D we assumed that the theoretical error of calculation of 
the muonium hyperfine splitting is about 70 Hz. This theoretical error is determined by the 
estimate of the still uncalculated terms which include single-logarithmic and nonlogarithmic 
radiative-recoil corrections of order a'^{Za){m/M)EF generated by the graphs containing 
besides the polarization loops also radiative photons, as well as the nonlogarithmic con- 
tributions of order {ZaY{m/M)Ep, a{Za)'^{m/M)Ep, and some other corrections (see a 
more detailed analysis in |^,|T^). Calculation of all these contributions and reduction of the 
theoretical uncertainty of the hyperfine splitting in muonium below 10 Hz is the current 
task of the theory. As the next step towards this goal we hope to present soon the results of 
calculation of the single-logarithmic and nonlogarithmic radiative-recoil corrections of order 
a'^{Za){m/M)Ep generated by the graphs containing besides the polarization loops also 
radiative photons. 
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APPENDIX A: AUXILIARY INTEGRALS 

All contributions to hyperfine splitting in the main body of this paper are written with 
the help of the function 

f{^,k)^mk) + U',ik), (Ai) 

where the standard auxiliary functions $„(A;) where introduced in ||10|| ) 



Hik) = Wi^,) - ^, (A2) 



$^(fc) = -^,W{^,) + \, (A3) 



and 



W^(U = Jl + ^-l> eM = /^'fc'- (A4) 

In terms of these functions all high-momentum contributions to hyperfine splitting may be 
represented as linear combinations of the standard integrals 

Vin.n = l j^i^r^krKik), (A5) 



where I = 1, m = 0,1,2 and n = 0, 1. Calculation of these integrals was described in 10 



and we present here only the results for the two integrals which were not calculated in |]TU 



Fi2o = I In' (2/i) - 2 In^ (2/i) + (y + 4) In (2/i) + C(3) - y - 4 + ^ InV, (A6) 



Vui = - 3 In' (2^) - - In^ (2^) - - + - j In (2/x) - -((3) " J^ " 4 " 3 ^^ ^- (^7) 



In Section |rV| we encountered the integral 



f dk i„ k[{viT¥ -k)-- (tv/in^ - p - 2)] / <i^- ;;,-;'/j) (A8) 

7r2 209 



18 432' 



which may be calculated by changing the integration variable 

2 

^ ~ k + ^A + k'^' 
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(A9) 



A number of integrals with the function R{v) (see eq.(p7D) used in Section are collected 
below 



dv 
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^t;^ln(l -t;^) = 2C(3)ln2 + ^C(3) + ^Li4^^ 
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